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SUMMARY 


An extension of the Tridiagonal Reduction (FEER) method in Level 17 of 
NASTRAN for complex eigenvalue analysis is described. As in the case of real 
eigenvalue analysis, the eigensolutions closest to a selected point in the 
eigenspectrum are extracted from a reduced, symmetric, tridiagonal eigenmatrix 
whose order is much lower than that of the full-size problem. The reduction 
process is effected automatically, and thus avoids the arbitrary lumping of 
masses and other physical quantities at selected grid points. The statement 
of the algebraic eigenvalue problem admits mass, damping and stiffness 
matrices which are unrestricted in character, i.e., they may be real, symmetric 
or unsymraetric, singular or nonsingular. 

The basic concepts underlying the method are summarized and special 
features, such as the estimation of errors and default modes of operation are 
discussed. In addition, the new user-information and error messages, and 
optional diagnostic output relating to the complex Tridiagonal Reduction method, 
are presented . 

Some numerical results and initial experiences relating to usage in the 
NASTRAN environment are provided, including comparisons with other existing 
NASTRAN methods for complex eigenvalue extraction. 


INTRODUCTION 


The complex Tridiagonal Reduction method is an extension of the FEER 
algorithm (^ast Eigenvalue Extraction _Routine) for real eigenvalue analysis to 
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complex, algebraic elgenproblem formulations. A specified number of eigen- 
values lying closest to a selected point in the complex plane are sought, as 
well as the associated eigenvectors. As in the case of real eigenvalue 
analysis (ref. 1), these eigensolutions are extracted from a synunetric, tri- 
diagonal eigenmatrix whose order is much lower than that of the full-size 
problem. In fact, the size of this canonical, reduced matrix is of the same 
order of magnitude as the number of desired roots, even if the discretized 
system model possesses thousands of degrees of freedom. The reduction process 
is carried out via an automatic algorithm requiring a finite number of steps. 
Thus, a basic weakness of methods requiring the lumping of masses and other 
physical quantities at arbitrarily selected degrees of freedom (refs. 2-4) is 
avoided in reducing the problem size. 

With regard to computational speed, the complex Tridiagonal Reduction 
method is somewhat slower than the Hessenberg method (refs. 5 and 6) for 
small problems (on the order of one hundred or less degrees of freedom) , if 
all the existing eigensolutions are to be calculated. However, it becomes 
more efficient than the Hessenberg method when the number of requested eigen- 
solutions is much less than the full problem size. Moreover, for much larger 
problems, the central memory requirement of the Hessenberg method exceeds the 
capabilities of most large computers, so that it becomes unavailable as a 
solution option. This limitation does not exist in the case of the Tri- 
diagonal Reduction Method . 

The complex Tridiagonal Reduction method employs a single initial shift 
point, and hence only one matrix decomposition is required for each neighbor- 
hood chosen in the complex plane. It therefore is more efficient than the 
complex Inverse Power method, which typically performs many shifts and de- 
compositions for each region selected. In addition, both the complex Inverse 
Power method and the complex Determinant method require that the user supply 
the length and width of rectangular regions in the complex plane, within which 
the eigenvalues are desired, as well as the number of estimated roots in each 
region. This can be burdensome to the user, who usually does not have 
enough advance insight to select these parameters intelligently. An improper 
choice (e.g., a strip too wide, or too small an estimate on the number of 
roots within the atrip) can lead to an inordinately large number of compu- 
tations or failure to extract any roots at all within the allotted machine 
time. These disadvantages are eliminated in the complex PEER method, where 
the user is only required to select points, closest to which a specified 
number of eigensolutions are desired. 

The theory and computational procedures for complex analysis depart from 
those of real analysis in the following major respects: 

1. Both left and right bl-orthogonal vectors must be created in the 
process of constructing the reduced tridiagonal matrix. 

2 . The reduced tridiagonal matrix, while symmetric in form, is, in 
general, complex rather than real, 

3. The calculated theoretical error.? in the cqmputed eigenvalues are 
estimates rather than upper bounds. 
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Eigensolutlons closest to one or more specified points (shift points) 
in the complex plane are found. All eigensolutlons obtained for 
previous shift points are swept out of the problem to prevent their 
regeneration when dealing with the current shift point. 

AN OVERVIEW OF THE CO^Q>LEX TRIDIAGONAL REDUCTION METHOD 

A detailed development of the analytical and computational procedures, 
including progranuning aspects and flow charts can be found in the NASTRAN 
Level 17 Theoretical and Programmer's Manuals. The following is a summary of 
the basic features of the complex Tridiagonal Reduction Scheme. 

The general complex eigenvalue problem is stated in the fonn 

[Mp“ + Bp + K]{u) = 0 , (1) 

where [M] , [B], and [R] may be real, complex, symmetric or unsymmetric, 
singular or non-singular. A specified number of eigenvalues, p , lying 
closest to a specified point, , (called a shift point) in the complex 

plane are to be found, as v»ell as the associated eigenvectors {u} . The 
eigenvalues may include multiplicities. By a suitable transformations the 
above can be expressed in the standard inverse form, 

(A]{x} = A{x} , (2) 

where [A] is double the size of the stiffness, mass and damping matrices, and 

A "1^ • 

* o 

In the special case where (B] is null (e.g,, no damping), the double-size 
eigenvalue problem can be avoided by considering the mathematical eigenvalues 
to be p- and defining 

A = (4) 

p -X~ 

in equation (2) . 

Since the eigenniatrix, [A] , is, in general, unsymmetric, the eigenvectors, 
{x} , are orthogonal to the eigenvectors, {x}, of the transpose eigenproblem 

[A]*^{x} = A{x} , (5) 

so that for A. A, , 

{xj}^'{x^} = 0; • (6) 

The above relationship is a biorthogonality condition and the associated 
eigenvectors, {x , } and {x } , are called right and left eigenvectors, 
respectively . ^ 

A reduction of the order of the eigenvalue problem, equation (2), is 
effected through the transformation 
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{x} = [V] {y} 
nxl nxm jnxl 

and 

{x} “ tv] {y} 
nxl nxm mxl 

where {x} and {x} are approximations of 
n Is the order of the unreduced problem, and m 4 n . The above trans- 
formation matrices are chosen to be blorthonormal, so that 

[V]^[V] = [I] . (8) 

From equations (2), (7), and (8), It is seen that 

[H]{y} « K{y} , (9) 

where 

[H] = [V]^[A][V] , (10) 

mxm 

and A is an approximation of the eigenvalue, A , 

Thus, equation (9) is an mth order eigenvalue problem, where m < n . 
The value of m is established according to the criteria given later. 

As in the case of real eigenvalue analysis (ref. 1), the Lanczos 
algorithm is used to construct the transformation matrices vector by 
vector, i.e.. 


(7a) 


’ (7b) 

{x} and {x} , respectively, 


[V] = [{v^}, {v^} , (11a) 

nxm 

[V] = [{v^}, {V2},....,{v^}] , (lib) 

nxm 

such that the reduced mxm matrix, [K], is tridiagonal and its liigenvalues 
accurately approximate the roots of equation (2) having the largest 
magnitude (or, equivalently, the physical roots, p , closest to the specified 
point of interest, Aq , in the complex plane) , Using symmetry arguments 
similar to those employed for real eigenvalue analysis, it can be shown 
that the transformed, reduced eigenmatrix in equation (9) is tridlagonal and 
symmetric , having the form. 
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nn « [v]'^'[A][v] - 



( 12 ) 


L ””™u 

The matrix coefficients are theoretically given by the simplified recurrence 
formulas 










1+1 




i “ l,m 


(13u) 


oaiOlNAL PACrlS 
POOR QUALPl 


i = 1, m-1 


(13b) 


where the sequence is initialized by choosing random, biorthonormal starting 
vectors for setting dj^-0;{v^]={v^}= {o} * 

The fin£il off-diagonal term, '^iu 4 -X * given by equations (13) is used in 
establishing error estimates for the computed eigenvalues, as described below. 
In addition, the above algorithm is modified in the computational scheme as 
follows: 


1. Each pair of vectors (v, .}, calculated in equations 
(13b), is reorthogonalizeJ to all previously computed pairs, 
before re-entering equations (X3a) . 


2. The size, m , of the reduced problem is a function of the number of 
accurate eigenvalues requested by the user and is limited to the num- 
ber of finite physical eigenvalues available. 
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The eigenvalues, K , and eigenvectors, {y}, of equation (9) are ex- 
tracted using the Q-R Iteration algorithm and eigenvector computational 
scheme described in connection with the Upper Hessenberg method in NASTRAN 
(ref. 5). They are then converted to physical form. 


CRITERIA FOR THE SIZE OF THE 
REDUCED EIGENVALUE PROBLEM 


The maximum number of finite eigensolutions, including any existing 
rigid body modes, is equal to the rank, r , of the eigenmatrix, [A], in 
equation (2). Thus, for example, massless degrees of freedom, appearing as 
zero diagonal terms in the [M] matrix, will result in singularities (rank 
reduction), which imply infinite physical eigenvalues. These spurious roots 
are swept out of the problem in the complex FEER process, with a consequent 
reduction in the available eigensolutions. 


A further consideration in limiting the maximum problem size is that the 
user has the option of requesting eigensolutions in the neighborhood of 
several shift points (Xoi,Xo2> • • •) in the complex plane. In the Tri- 
diagonal Reduction method, all eigensolutions, f , obtained for previous 
shift points are swept out of the problem to prevent their re-genera ion when 
dealing with the current shift point. This implies that the maximum possible 
size, m , of the reduced problem is further limited to 


m = r 
max 


- f 


(14) 


On the basis of numerical experiments, similar to those cited in 
reference 1 for real eigenvalue analysi"*, it has been found that when 
m << m , a first grouping of more than m/2 computed eigenvalues closest 
to the s^iift point are in accurate agreement with the corresponding number 
of exact eigenvalues, provided that 7 < m ^ . The remaining reduced- 

system roots are spread across the remaining exact elgenspectrura. To enhance 
the accuracy of the associated eigenvectors, the minimum problem size is 
further increased to twelve, again assuming that n* Wmax * 


Thus, if the user requests a total of q eigenvalues closest to a 
specified point in the complex plane, the order of the reduced problem is 
initial 1"' set to 

f =min[(2q+10),(2n-f)]; •B]?t[o] , (15a) 

\ 

( =min[ (2q+10) , (n-f) ] ; [B]=[0] . (15b) 


Although the total number of eigensolutions requested should not exceed 
%iax > there is usually no simple way to discern this upper limit in complex 
eigenvalue problems. However, the reorthogonalization tests are designed to 
automatically establish this upper limit. If the latter tests fail for some 
vector pair ’ this is an indication that a null vector has 
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been generated because mmax linearly independent vectors have already been 
obtained. The recurrence algorithm, equations (13), is then terminated and 
the order of the eigenproblem is further reduced to m = i . 


ERROR ESTIMATES FOR THE 
COMPUTED EIGENVALUES 


Following a development similar to that of reference 1 for real eigen- 
value analysis, it can be shown that 






( 16 ) 


The above shows that the absolute value of the difference between the com- 
puted and true eigenvalue magnitudes is proportional to the magnitude of 
dm+1 (which is the next off-diagonal term that would be generated had the 
reduced tridiagonal matrix, [H] , beer: increased from order m to order 
m+1 ) and y„jj^ , (which is the last term in the reduced-system eigenvector 
associated with A ) . 


Converting equation (16) to physical eigenvalue form, using equations 
(3) and (4), yields, 


^i 


(W 


- 1 


^ to] , 
l^il 


(17a) 




i 



- 1 


lAj 


[B] = [0] . 


(17b) 


The use of the above error estimates as criteria for acceptable eigensolutions 
is as follows: 


(a) If the physical eigenvalue, , corresponds to a zero root (e.g,, a 

rigid body mode) , the above computational scheme is invalid and there- 
fore bypassed* Denoting t as the number of digits carried in the 
computations, a zero root is assumed to occur whenever 


where 



RMS 


< 1 


.-t/3 


9 


RMS -i [1J2| 


1 - 2 , 


+ 




and is denoted by setting the error 


to zero . 


(18) 

(19) 
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(b) The eigenvalues are listed in order of increasing distance from the shift 
point, Xq , to determine whether their associated estimated errors, , 
meet an acceptable relative error tolerance set by the user on the EIGC 
bulk data card (the default value is 0.10/n, where n is the order of 
the stiffness matrix) . The first eigenvalue not meeting the tolerance 
test, as well as all subsequent eigenvalues further removed from the 
center of interest, are considered to lack sufficient accuracy and are 
therefore discarded. 

(c) Acceptable eigenvalues obtained in the above manner are reordered 
according to the magnitude of the imaginary part, with positive values 
considered as a group ahead of all negative? values. 


NASTRAN USER'S INSTRUCTIONS 


Figure 1 shows modifications of the EIGC card in the NASTRAN bulk data 
deck which accommodate user implementation of the Tridiagonal Reduction method 
for complex eigenvalue analysis. The modifications consist of additions to 
the standard user instructions and are underscored for ease in identification. 

When the complex Tridiagonal Reduction method is invoked, the E 
parameter on this card represents the maximum allowable value of the com- 
puted absolute relative error in a physical eigenvalue. If this value is 
exceeded, the associated eigensolution is not accepted for further processing 
by NASTRAN. A detailed list of the maximum relative errors computed by com- 
plex PEER can be obtained by requesting DIAG 12 in the NASTRAN Executive 
Control Deck. 


USER MESSAGES AND OPTIONAL DIAGNOSTICS 


Functional Module User Messages 

The following is a description of the NASTRAN user messages which may 
be generated by NASTRAN during the execution of the Complex Tridiagonal 
Reduction method and which are unique to this method . Explanatory infor- 
mation is provided following the text of each message and, in the case of 
a fatal message, corrective action is indicated. Refer to the NASTRAN 
Users' Manual, Section 6 for a complete listing of other system and user 
messages . 

Fatal messages cause the termination of the execution following the 
printing of the message text. These messages will always appear at the 
end of the NASTRAN output . Warning and information messages will appear at 
various places in the output stream. Such messages convey only warnings or 
information to the user. Consequently, the execution continues in a normal 
manner following the printing of the message text . 
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Message List 


3149 *** USER WARNING MESSAGE 3149, USER SPECIFIED NEIGHE0Rli00D CENTERED 

AT 0RIGIN N0T ALL0\ffiD, CENi'ER SHIFTED T0 THE RIGHT .001. 

Poin^. of interest in the complex plane (oiai>^ai) » closest to 
which the eigenvalues will be computed, was input as (0,0, 0.0) 
on an EIGC bulk data continuation card. Since this is an In- 
admissible choice, the point automatically used was (.001, 0.0). 

3150 *** USER WARNING MESSAGE 3150, DESIRED NUMBER 0F EIGENVALUES 

*** INVALID. SET = 1. 

Number of accurate roots desired Ndl , was omitted, input as aero 
or negative on an EIGC bulk data continuation card. The number 
automatically used was 1. 

3151 *** USER WARNING MESSAGE 3151, DYNAMIC MATRIX IS SINGULAR 

(0CCURRENCE ****) IN NEIGHB0RH00D CENTERED AT **** **** 

Point of interest in the complex plane (cx^isWai), closest to 
which the eigenvalues will be computed, was input too close to 
an eigenvalue on an EIGC bulk data continuation card. The point 
is automatically shifted by adding .02 to both the real and 
imaginary parts. If the dynamic matrix is still singular, the 
next neighborhood, if any, is searched. 

3152 *** user' INE 0RMA ■ i)N MESSAGE 3152, SUBR0UTINE ALLMAT 0UTPUT 

EIGENVALUE IS NULL. 

Wien an eigenvalue output from subroutine ALLMAT is exactly zero, 
the formula for computing the associated theoretical error test 
fails. The magnitude of the eigenvalue is considered to be 10~10 
for use in that formula. 

3153 *** USER WARNING MSSAGE 3153, ATTEMPT T0 N0RMALIZE NULL 

VECT0R IN SUBR0UTINE CFEER4. N0 ACTION TAKEN. 

An eigenvector output from subroutine ALLMAT is a zero-vector. 

315A USER WARNING MESSAGE 3154, SIZE 0F REDUCED PR0BLEM 

DECREMENTED 0NCE (N0W ****) DUE T0 NULL ERR0R ELEMENT. 

If subroutine GFEER4 receives a reduced tridiagonal matrix having 
error element duj^-i exactly (0,0), it is impossible to compute 
meaningful theoretical error estimates for any of the eigenvalues. 
The size of the reduced problem is reduced by one, so that d^ 
becomes the new error element. 
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3155 *** USER WARNING MESSAGE 3155, REDUCED PR0BLEM HAS VANISHED. 

N0 R00TS T0UND. 

If decrementinfe the size of the reduced problem (see message 
315^) causes the. size to become zero, the program continues to 
the next neighborhood, if any. 

3156 *** USER WARNING MESSAGE 3156, SIZE 0F REDUCED PR0BLEM 

REST0RED T0 **** BECAUSE NEXT ERR0R ELEMENT WAS ALS0 
NULL. ERH0R ELEMENT SET = **** **** 

This message follows message 3154. If d^ is also exactly zero 
(in addition to dj^l being exactly zero), then the original re- 
duced problem size is restored and is set to (e,0) where 

£ = E/lOO and E is the error tolerance on acceptable eigenvalues 
input on the EIGC bulk data card. 

3157 ***■ USER WARNING MESSAGE 3157, PEER PR0CESS MAY HAVE 

C;u.CULATED FEWER ACCURATE M0DES **** THAN REQUESTED 
IN THE NEIGHB0RH00D 0F **** **** 

The desired number of eigenvalues specified on the EIGC bulk 
data continuation card exceeds the additional number that can 
be calculated by the Complex Tridiagonal Reduction (Complex 
FEER) method in the current neighborhood . 

3158 *** USER WARNING MESSAGE 3158, N0 ADDITI0NAL M0DES CAN BE 

F0UND BY FEER IN THE NEIGHB0RH00D 0F **** **** 

An initial pseudo-random vector cannot be made orthogonal to 
the existing set of orthogonal vectors (which come from Restart 
and from all prior-neighborhood sets of eigensolutions) . 

3159 *** USER INF0RMATX0N MESSAGE 3159, ALL S0LUTI0NS HAVE BEEN 

F0UND. 

The FEER method has salved the entire problem. Any additional 
neighborhoods (as specified by the presence of EIGC bulk data 
continuation cards) are ignored. 

3160 *** USER INF0RMATION MESSAGE 3160, MINIMUM 0PEN C0RE N0T 

USED BY FEER **** W0RDS {***^ K BYTES) . 

This message indicatet; the amount of open core, in both bytes 
and words, not used by FEER. 

3161 *** USER WARNING MESSAGE 3161, DESIRED NUMBER 0F EIGENS0LU- 

TI0NS **** F0R NEIGHB0BH00D **** 0F **** CENTERED AT 

AAA* exceeds the EXISTING NUMBER ****, ALL EIGENS0LU- 
TI0NS WILL BE S0UGHT,, 


428 



The desired nvuuber of eigenvalues specified on the EIGC bulk data 
continuation card exceeds the size of the eigenmatrix, which is 
the maximum possible number of existing eigenvalues. 

3162 A** USER WARMING MESSAGE 3162, ATTEMPT T0 NORMALIZE NULL 

VECTOR. N0 ACTI0N TAKEN. 

The general vector normalization routine (CFN0R1 or CPK0R2) has 
a zero-vector input to it. 

3163 *** USER WARNING MESSAGE 3163, ALL **** S0LUTI0NS HAVE FAILED 

ACCURACY TEST. M0 R00TS F0UND. 

The number of eigensolutions passing the relative error test is 
zero. The maximum allowable error for the relative error test 
is specified in field 7 of the EIGC bulk data card, A detailed 
list of the computed error bounds could have been obtained by 
requesting DIAG 12 in the Executive Control Deck. 

316-!i *** USER INF0RMATI0N MESSAGE 3164, ALL **** S0LUTIONS ARE 
ACCEPTABLE . 

All the eigensolutions obtained in the reduced problem correspond- 
ing to the point of interest pass the relative error test. The 
maximum allowable error for the relative error test Is specified 
in field 7 of the EIGC bulk data card. A detailed list of the 
computed error estimates could have been obtained by requesting 
DIAG 12 in the Executive Control Deck. 

3165 *** USER INF0RMAT10N MESSAGE 3165, **** S0LUTI0NS HAVE BEEN 

ACCEPTED AND **** S0LUTI0NS HAVE BEEN REJECTED. 

Some eigensolutions pas.sed the relative error test and some 
did not. 

3166 *** USER INF0RMATI0N MESSAGE 3166, **** l't0RE ACCURATE EIGEN- 

S0LUTI0NS THAN THE **** REQUESTED HAVE BEEN F0UND F0R 
NE1GHB0RH00D **** 0F **** CENTERED AT **** AAAA. USE DIAG 
12 T0 DETERMINE ERR0R ESTIMATES. 

The number of eigensolutions passing the relative error test is 
greater than the number requested on the corresponding EIGC bulk 
data continuation card. The maximum allowable error for the 
relative error test is specified in field 7 of the EIGC bulk 
data card. A detailed list of the computed error estimates 
could have been obtained by requesting DIAG 12 in the Executive 
Control Deck. 
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The Eii^envalue Summary Table 

The following summary of the eigenvalue analysis performed, using the 
complex Tridiagonal Reduction (FEER) method, is automatically printed: 

1. Number of eigenvalues extracted, 

2. Number of starting points used. 

This corresponds to the total number of random starting and restart 
vectors used by the complex FEER process for all neighborhoods, 

3. Number of starting point moves. 

Not used in FEER (set equal to zero), 

4. Number of triangular decompositions. 

Always equal to the number of poin^. g of interest (neighborhoods) 
in the complex plane processed by FEER, since ordinarily only 
one triangular decomposition is required by FEER for each point of 
interest, unless the dynamic matrix is singular at a given point 
of interest, in which case an additional decomposition is required 
(obtained by moving the point of interest slightly) . 

5. Total number of vector iterations. 

The total number of reorthogonalizations of all the trial vectors 
employed . 

6. Reason for termination, 

(.0) All, or more solucions than the number requested by the user, 
have been determined (normal termination) . 

(1) All neighborhoods have been processed, but FEER has not obtained 
the desired number of roots in each neighborhood, possibly be- 
cause they have already been found in other neighborhoods, 

(2) Abnormal termination - either no roots found or none pass the 
FEER error test* 


Optional Diagnostic Output 

The user can obtain special detailed information relating to the 
generation of the reduced problem size, the elements of the reduced tri- 
diagonal matrix, vector reorthogonalization iterations, computed error 
estimates, order of eigenvalue extraction, and distance of extracted 
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elgenvalaci from the ceutei' of interest by requesting DIAG 12 in the NASTRAN 
executive control deck. 

The meaning of this infoi^niation is explained below in the order in which 
it appears in the DIAG 12 output. 

(FAST EIGENVALUE EXTRACTI0N R0UTINE) 

This header is always prin'^.-d first. 
aaaA'SXNGLE precision W0RDS 0F 0PEN C0RE N0T USED (SUBR0UT1NE XXXX) 

•kicii-k - Open core not used by subroutine XXXX, in single-precision words, 

XXXX - Either CFGNTL, CFEER3, or CFEER4 . This message appears 
three times* 

CFCNTL ACCURACY GRITERI0N * (INPUT VALUE**) 

A - Accuracy criterion, used for rejecting eigensolutions (expressed 
as a percentage) . 

A A Value ot accuracy criterion input by the user on the EIGC bulk 

data card. 

CFCNTL NEICHB0Rli00D * CENTER = ** kk N0.DES.RTS. = *** N0NSYM = A*** 

A - Positive integer indicating which neighborhood, or center of 
interest, is currently being processed. 

kk kk - Center of interest in the complex plane. 

AAA - Number of desired roots for the current neighborhood, input by 
the user on the corresponding EIGC bulk data continuation card. 

kkitiK - Indicator which, when nonzero, forces the program to consider 
the matrices as non-'synu«edric, even though they may actually 
be symmetric. This is input by the user in field 7 of each 
EIGC bulk data continuation card. This input was used during 
program checkout of the complex FEER process, and it should 
have no affect on the solution* However, the user should leave 
field 7 blank on each EIGC continuation card. 

RE0RT110G0NALIZATI0N ITE1L\TI0N * TARGET VALUE == *A 
ERRORS = AAA AAA AAA AAA 

A - The reorthogonalization iteration-number. This message will 
appear many times, as the FEER process **cleans up" each trial 
vector by forcing it to be as orthogonal as possible to the set 
of vectors already computed. 

A A _ Convergence tolerance, such that the errors must be smaller than 

this value* In order to avoid taking square roots, the tolerance 
and errors are all squared. 

AAA AAA AAA AAA Four reorthogotializa tion errors, the first two of 

which correspond to the orthogonality of the current right and 
left handed trial vectors, respectively, with respect to all 


431 



previously computed vectors In the current neighborhood, and the 
latter two of which correspond to the orthogonality of the same 
vector's with respect to all eigenvectors previously computed 
(restart and prior neighborhoods) . 

REDUCED TRIDIAGCiNAL MATRIX ELEtlEWTS R(^iW * 

0FF diagonal « ** 

DIAGONAL ~ 

* - The row number of the (reduced) tridiagonal matrix. 

A'lV :k:k - Value of the of f-'diagonal element for that row. 

-^ick « Value of the diagonal element for that row. 

Following the printing of several lines containing reorthogonalization 
information and reduced tridiagonal matrix elements, when the PEER process 
has finished its computations for the current point of interest, the header 
(see above) is printed once again, follox^ed by a table which summarizes all 
the eigensolutions found by PEER. This table has seven columns, as follows: 

(1) Solution number. This is simply a positive integer 1,2,3,*.. 

(2) Order of extraction. These numbers indicate the order in which 
the tridiagonal matrix v;as constructed, 

(3) Distance from center. This is the distance from the extracted 
eigenvalue to the neighborhood center (which is printed above 
the table) in the complex plane. The tabular values are sorted 
according to increasing distance from the center. 

(4) Real part of the extracted eigenvalue* 

(5) Imaginary part of the extracted eigenvalue. 

(6) Theoretical error estimate. This value must be smaller than the 
Accuracy Criterion (see above) for the eigensolution to be 
acceptable . 

(7) Status. A single word, ^'accept^' or "reject*^, to indicate the 
result of the accuracy test. A minus sign (-) is added to 
‘'reject*' so that the eye can more rapidly distinguish between 
the two words. 

Finally, this table is printed a second time, but with the rejected 
eigensolutions deleted* 

For very small problems, there is a Very Detailed Printout (VDP) option. 
This option was originally used to debug the complex PEER logic, and is no 
longer required. DIAG 12 must be specified in the executive control deck to 
invoke the VDP option, and furthermore, field 6 of a given EIGC bulk data 
continuation card must have a (floating point) value equal to or greater than 
the size of the stiffness matrix* Thus, field 6 of each EIGC continuation 
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card should ordinarily be left blank, so that the VDP option will be 
suppressed. The actual printed output of this option consists of all 
vectors for each step of the coiiiplex PEER process, which is too extensive and 
detailed for normal user purposes, 

NUMERICAL RESULTS AND COMPARISON WITH 
EXISTING NASTRAN EIGENVALUE METHODS 

Test Problem 

During the developmental stages of the complex Tridiagonal Reduction 
method, a simple three degree of freedom model consisting of the rod and 
dashpot system shown in figure 2, was employed for experimental and check- 
out purposes. The eigenvalues of this system consist of one complex- 
conjugate pair, one pure imaginary and three zero roots,, the latter 
corresponding to rigid-body modes. 

A comparison of the performance of the complex Tridiagonal Reduction 
method, the complex Determinant method and the complex Inverse Power method, 
is summarized for this example in table I, In run 1, all six eigensolutions 
were successfully found using complex PEER and a single shift point in the 
second quadrant of the complex plane; only one triangular decomposition of 
the dynamic matrix was required. In run 2, using the complex D etem\inant 
method, 29 triangular decompositions were performed and only two, non-zero, 
eigenvalues were found, in spite of the fact that muHtiple search regions 
were used, including one region encompassing the origin. The reason given 
for termination in the eigenvalue sumn\ary table was that ''all predictions 
for eigenvalues are outside the regions specified'*, even though this was 
not the case. 

Runs 3-5 were with the complex Inverse Power method, which did not fare 
too well in this exercise. A large number of iterations were perfomued, re- 
sulting in relatively high CPU and I/O times. In addition, the first run 
yielded only two roots and the remaining two, none at a' I, even though 
the search regions selected, while differing from run to run, encompassed 
the knoiim eigenvalues . 

All the above runs, and those reported below, were performed on the 
NASA/GSFC IBM 360/95 computer. 


Tidal Frequencies and Modes in Closed, Shallow Basins 

CO! current with the complex PEER development, an independent study was 
conducted (ref. 7) to estimate the tidal frequencies and mode shapes in two 
of the Great Lakes, namely, Lake Erie and Lake Superior. It was decided to 
use a finite-element displacement fomiulation and the complex eigenvalue 
extraction capabilities in NASTRAN for this purpose. The mathematical model- 
ing was based on the following assumptions: 

(a) Negligible convective accelerations and other nonlinear effects. 
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(b) Inviscid, incompressible flow, 

(c) Negligible vertical velocity compared to lateral velocities. 

(d) The Boussinesq hypoUheaio of hydrostatic pressure variations, 

(e) Sufficient shallowness to permit integration-averaging through 
the depth. 


As a consequence of the above simplifications, the Navier-Stokes equations 
reduce to the two-dimensional form, 



|~ + S 1^ " fv = 0 

(20a) 


8v , 3C j. « 

+ 8 5^ + £u . 0 

(20b) 


If + (hu) + (hv) = 0 , 

(20c) 

where 



u,v 

? 

h 

f 

g 

= integrated average velocities in the x and y (horizontal) 
directions 

= local water height, measured from the mean surface 
^ local depth of water, measured from the mean surface 
= the Coriolis parameter, 2f2sin(|) , with (j) the latitude, and 
the earth’s rotational rate 
= acceleration due to gravity. 


Using the Galerkin method, a finite elenient representation of the above 
equations was developed, in which the nodal variables are u , v and C • 

A data-generator code was then written, which generated DMIG card images for 
use as NAS TRAN input * 

The two lakes are shown in figures 3 and 4 and their finite-element 
meshes are given in figures 5 and 6- The mesh of Lake Erie contains 81 nodes 
and 204 unconstrained degrees of freedom, while the representation of Lake 
Superior involves 124 nodes and 299 unconstrained degrees of freedom. These 
selections were based on numerical convergence studies with successively 
finer meshes, and represent the fineness needed to obtain two or three 
accurate modes . 

In the early stages of the study, the complex Inverse Power method was 
used, but had to be abandoned because of inconsistencies in the results; 
the roots obtained seemed to depend on the search region selected and false 
roots were almost always calculated very close to the starting point in the 
region, A tightening of the convergence criterion on the EIGC bulk data 
card was attempted, but this did not resolve the difficulties. 
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Some limited success was achieved using the complex Doi.enni: * .C method, 
but here again, difficulties were encountered* Unless the search region 
could be rather carefully and precisely defined, the chances of finding roots 
was somewhat remote* The method is apparently quite sensitive in its search 
pattern, and does not home-in on a root if there is an extensive search area 
to work on* 

At some point approximately mid-way in the study the complex PEER 
capability was completed and integrated into NASTRAN, whereupon it was 
applied to the tidal mode problem for the lakes* The results obtained were 
consistently good. To check their accuracy, changes were made in the shift 
points, mathematical scaling and the number of requested eigensolutions 
(i.e*, the truncated size of the probiem) * These variations had only a 
negligible effect on the calculated frequencies and mode shapes. 

Some timing results for Lake Erie, using complex PEER and the Determinant 
method are given in table II. As noted above, tlie Determinant method was 
viable only if the search regions were made very small, implying that the 
locations of the roots were rather well kno™ in advance* 

The Upper Hessenberg method could not be used for these models, since 
the maximum region available on the IBM 360/95 was 900K., which would only 
permit a 43 node mesh with this metliod* 


CONCLUDING REMAPJCS 


Initial experiences with the complex Tridiagciuil Rocluction (PEER) 
method indicate that it is very effective in extracting ‘.ny desired number 
of accurate complex eigensolutions in the neighborhood of a selected shift 
point on the complex plane. The method automaticalJ y computes complex 
roots at increasing distances from the selected poiv^t until the requisite 
number, specified by the user, is obtained. In this respect, a disad- 
vantage of the complex Determinant and complex Inverse Power methods, 
namely, a very careful delineation of search regions, is eliminated* In the 
case of multiple shift points, it has been found that complex PEER success- 
fully sweeps-out eigensolutions obtained for previous shift points and 
prevents their regeneration when dealing with the current shift point. 

Since the mathematical properties and characteristics of complex 
eigenvalue problems are very broad and varied, it should be recognized that 
the results reported herein with regard to computational efficiency and 
timing, are only indicative of a small class of problems* A fuller assess- 
ment of the capabilities of this new method can only be obtained following 
extensive application experiences within the user community. 
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TABLE I. COMPARISON OF EIGENVALUE METHODS 
FOR ROD-AND-DASHFOT TEST PROBLEM 


Run 

Mechod 

I339I 

mm 

Region 

CK) 

Comments 

1 

Complex 

FEER 

i . 

0.151 

2.262 

300 

Eigensolutions Found: All 6 

No. of Triangular DGcompositions: 1 
Search Regions: One point in 2nd 

quadrant of complex plane 
Reason for Teniiination; All solutions 
found , 

2 

1 

Complex 

Determinant 

0.174 

2.858 

320 

Eigensolutions Found: 2 (3 rigid body 
modes at origin missed) 

No. of Triangular Decompositions: 29 
Search Regions: 3, around knowi roots, 
including origin 

Reason Given for Termination: All 
predictions are outside regions 
specified 

3 

Complex 1 
Inverse 
Power 

0.271 

8.997 

320 

Eigensolutions Found: 2 (3 rigid 
body modes missed) 

No. of Triangular Decompositions: 1 
in last search region 
Search Regions: 3, around known 
roots, including origin 
Reason Given for Termination: Number 
of desired roots have been found 

h 

i 

Complex 

Inverse 

Power 

! 

0.238 

i 

4.100 

1 

320 

1 

Eigensolutions Found : None 

No. of Triangular Decompositions: 4 
in last search region 
Search Regions: 1, encompassing all 
5 existing roots on and above real 
axis 

Reason Gxven for Termination: Four 
starting point moves while track- 
ing a single root 

5 

Complex 

Inverse 

Power 

0,149 

1,566 

320 

1 

Eigensolutions Found: None 

No- of Triangular Decompositions: 2 
Search Regions; Around origin, in 
attempt to find rigid-body modes 
Reason Given for Termination; Two^ 
successive singularities found 
while performing triangular 
decomposition 
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TABLE II. TIMING RESULTS FOR LAKE ERIE TIDAL MODE 
CALCULATIONS - 81 NODE MODEL 


Run 

Method 


Bn 

■H 

Comments 

1 

1 

Complex 

PEER 

2.164 

! 

5.886 

700 

6 accurate modes requested 
4 obtained 

2 

2.425 

8.593 

350 

15 accurate modes requested 
12 obtained 

3 


7.576 

27.751 

560 ! 

50 accurate modes requested 
54 obtained 


Complex 

6.572 

1 

8,618 

700 

3 modes obtained; insufficient time 
for more 

5 

Determinant 

5.315* 

7.489 

700 

6 modes obtained; insufficient time 
for more 


*Less time used due to a more careful choice of the search region. 
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OHIGINAL PAGE IS 
OP POOR QUALITYl 


BULK DATA DECK 


Input Data Card EIGC Complex Eigenvalue Extraction Data 
Description ; Defines data needed to perform complex eigenvalue analysis 





Format and Example ; 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


SID 

METHOD 

N0RM 

G 

C 

E 1 



+abc 

1 EIGC 

1 14 

DET 

P0INT 

27 


1.-8 



ABC 








\l 

^dl 


+def 

+BC 

2.0 

5.6 

2.0 

-3.4 

2.0 

4 

4 

ZZZ] 

DEE 


-{-def 

“a2 

a 2 

%2 

^^2 

^2 

«e2 

«d2 



+EF 

-5.5 

1 

Ul 

5.6 

5.6 

1.5 

6 

3 




(etc • ) 


Figure 1. Modifications to the EICC bulk data card for the Tridiagonal 

Reduction Method. 
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Field 


Contents 


SID 

METHtiD 


M0SM 


G 

C 

E 


'“ar “ajM 


Set identification number (unique integer > 0) 

Method of complex eigenvalue extraction, one of the BCD values, 
"INV", "DET", "HESS” or " FEER" 

INV “ Inverse power method 

DET - Determinant method 

HESS - Upper Hessenberg method 

FEER - Tridiagonal Reduction Method 


Method for normalizing eigenvectors, one of the BCD values 
"MAX” or "P0INT" 

MAX - Normalize to a unit value for the real part and a 
zero value for the imaginary part, the component 
having the largest magnitude 

P0INT - Normalize to a unit value for the real part and a 

zero value for the imaginary part the component de- 
fined in fields 5 and 6 - defaults to "MAX" if the 
magnitude of the defined component is zero. 

Grid or scalar point identification number (Required if and 
only if N0RM«P0INT) (lnteger>0) 

Component number (Required if and only if N0RM="P0INT" and G 
is a geometric grid point) (0^^integer>6) 

Convergence criterion (optional) (Real^O . 0) 

F or method = "FEER", error-tolerance on acceptable eigenvalues 
in percent (default value is .IQ/n, where n is the order of 
the stiffness matrix) 


Two complex points defining a line in the complex plane (Real) 
For method = "FEER" , (a^^ j , ) is a point of Interest in the 

complex p lane, closest to which the eigenvalues are computed ; 
|g^. I + > Q * The point is ignored . 


a. 

J 


Width of region in complex plane (Real>0.0) 
Blank for method = "FEER" . 



Estimated number of roots in each region (lnteger>0) 
Ignored for method ~ "FEER" . 


Figure 1 . Continued 
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N 


dj 


Desired number of roots in each region (Default is 3Ngj) 
(lnteger>0) Desired number of accurate roots for method " 
"FEER” (Default is 1) . 


R'^marks ; 

1. Each continuation card defines a rectangular search region. For 

me thod = *'FEER", the card defines a circular search region, centered 
at (otaj 1 Waj) and of sufficient radius to encompass roots . 

Any number of regions may be usea and they may overlap. Roots in 
overlapping regions will not be extracted more than once. 

2. Complex eigenvalue extraction data seta must be selected in the Case 
Control Deck (CMETH0D=SID) to be used by NASTRAN. 

3. The units of a, w and S. are radians per unit time. 

4. At least one continuation card is required. 

5. For the determinant method with no damping matrix, complex conjugates 
of the roots found are not printed. 

6. See Section 10.4.4.5 of the Theoretical Manual for a discussion of 
convergence criteria, 

7. For the Upper Hessenberg method, controls t! .ber of vectors 

computed. Only one continuation card is considereu and the (a, w) 
pairs, along with the parameters and , are ignored, ln~ 

sufficient storage for HESS will cause the program to switch to INV. 

8 . The error tolerance, E , for the "FEER” method is with regard to 




- 1 


for [B] i= [0] and 


|p^-(a , , w 

'*^1 a,i* a.r ‘ _ ^ 

|p?-(a OJ .)^1 

' 1 a j aj ‘ 


for [B] = [0] 


where is a computed eigenvalue and an exact eigenvalue, 


Figure 1 . Concluded . 
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DASHPOT 



Figure 2, 


Test Problem — Rod and Dashpots with 3 Degrees of Freedom 





Figure 3^ Illustration of Lake Erie and the adjacent geography- 
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Figure 5. 81 Node Finite Element Model for Lake Erie, 
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Figure 6. 


124 uode finite element -model for Laite Superior. 


